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Abstract 

In this paper a system-oriented formalism of Quantum Information 
Processing is presented. Its form resembles that of standard signal 
processing, although further complexity is added in order to describe 
pure quantum-mechanical effects and operations. Examples of the ap- 
plication of the formalism to quantum time evolution and quantum 
measurement are given. 

1 Introduction 

The field of Quantum Information Processing (QIP) has recently received 
a renewed attention, particularly by the Physics and Computing commu- 
nities [[l], H|, H|, ||, [|. This discipline, however, dates back from the six- 
ties and seventies, when the advent of the laser generated a lot of interest 
among communication-theory researchers (see ||, |?J and references therein). 
This interest arose from the possibility of employing the laser as an elec- 
tromagnetic source in ultra-wideband communication systems. In order to 
establish a suitable model of the communication channel, and owing to the 
high frequency of the laser radiation, it was necessary to consider in the 
modelling the quantum-mechanical properties of radiation (quantum noise). 
These early research efforts constituted the birth of Quantum Communica- 
tion Theory. 

First explorations on the quantum nature of information are due to 
MacKay ||, although it was Stern || the first, to our knowledge, who mod- 
elled optical communication channels quantum-mechanically. Gordon JlO| ] 
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also recognized the influence of quantum noise in communication systems, 
and addressed the problem of selecting the appropriate receiver structure 
for an optical communication channel. Following the standard approach of 
Quantum Field Theory, She [11] studied the channel and the receiver of a 
communication system. Later, in what we consider was the most complete 
research effort towards the full quantum-mechanical characterization of the 
detection problem in an optical communication channel, Helstrom et al. || 
analyzed in detail how the laws of Quantum Mechanics (QM) affect the reli- 
ability of a communication system, and how to design quantum-mechanical 
receivers. These same issues were also addressed in a similar manner by Liu 
[l|, Personick g|], and Davies @, ||. 

During the last twenty-five years work has continued on the character- 



ization of quantum receivers and its optimum design JbJ, |l_y, |17|, £L8| [L9f| ; 
however, the full modelling of information-processing quantum-mechanical 
systems with a broader scope (not limited to optical communication chan- 
nels) has only recently been proposed (see, e.g. ]^0| and references therein), 
arguably fueled by the fertility in the field of Quantum Optics during the 
eighties. As a consequence, brand-new disciplines related to the quantum 
processing of information, such as Quantum Computation, Quantum Cryp- 
tography and Quantum Teleportation, have appeared ||]. Behind these new 
disciplines we find again QM, but unlike in the typical formulations of the 
sixties and seventies, now the processing (in communication as well as com- 
puting) is described making use of the theory of quantum open systems 

II- 

QM is not an easy theory to grasp, and Open-Systems QM is even more 
subtle, particularly to the Information-Processing community. With this 
paper we expect to contribute to shorten the gap that separates researchers 
with a Physics background from those more concentrated on classical Sig- 
nal Processing (SP). With this purpose in mind, we show that a general 
formalism of QIP similar to that of classical SP can be established. 

The paper is organized as follows. In Section § we review the formalism 
of classical SP. In Section |3| we define the concept of quantum system and 
discuss how to efficiently model the dynamics of such a system in general. 
In Section |] we derive the matrix representation of the input-output rela- 
tionship of a quantum system. Section || generalizes the proposed formalism 
to composite systems. This allows to describe the concept of entanglement. 
In Section ^ we apply the formalism to the description of two key quan- 
tum processes: Quantum Evolution and Quantum Measurement. Finally, 
in Section ^ we summarize the main results of the paper. 
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2 Classical signals and systems 



SP is concerned with the representation, transformation, and manipulation 
of signals and the information they contain. Generally, this information 
is about the state or behaviour of a physical system, although we shall 
regard signals just as containers of information. Mathematically, signals are 
represented as functions of one or more variables belonging to a vector signal 
space. 

Depending on the character of these variables, continuous or discrete, 
SP has been divided into Analog and Discrete Signal Processing (DSP), 
respectively. Progress in computer technology has made possible the easy 
implementation of the latter and its pervasiveness. Thus, in the following 
we shall restrict ourselves to DSP. 

Discrete signals are represented by sequences of numbers that we shall 
denote, following the standard formalism of one-dimensional DSP, by x[n] 
(n = 1, • • • , N). The actual processing of signals is performed by systems, 
that can be seen as entities that map an input signal x[n] to an output signal 
y[n]. A particularly important class of systems is formed by those that are 
linear and time- invariant. These two properties allow to write the system 
input-output relationship by means of the convolution operator: 

N 

y[n] = x[n] * h[n] = ^ x[k]h[n — ft], (1) 
k=l 

where the sequence h[n] is the impulse response of the system. The rela- 
tionship expressed by equation ([l]) can also be written in the transformed 
domain by means of the Discrete Fourier Transform (DFT): 

Y[k] = X[k]H[k], (2) 

where the capitalized sequences are the DFTs of the corresponding signals. 

Besides linearity and time-invariance, other restrictions of more funda- 
mental physical nature, such as stability and causality, can be imposed on 
the operation of systems. Further information regarding these and similar 
issues can be seen in [^] . 

3 Quantum signals and systems 

Just like a classical system, a quantum system is a mathematical entity 
that operates on the quantum signal at its input and generates an output 
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quantum signal. Quantum signals are often called states, and, following the 
standard QM formalism, are represented by vectors belonging to a Hilbert 
signal space £. In "bra-ket" or Dirac notation, we shall write the input and 
output signals of a system as \x) and \y), respectively. We shall assume 
in the following that the signal space has a finite dimension N. This is a 
requirement in order to describe discrete quantum signals. 

Quantum systems, as classical ones, can manipulate two types of sig- 
nals, pure and mixed (stochastic). Pure signals are deterministic; on the 
contrary, mixed signals result from the statistical mixture of pure ones. As 
an example, consider the case in which the input signal can be one of the set 
i = 1, • • • , M} with probability pi. In situations like this, the informa- 
tion about the input of the system can be better described by an hermitian 
operator called density operator (see ]^3| for a definition of density operators 
and their more relevant properties) that we write as: 

M 

Px = ^2pi\xi)(xi\. (3) 
1=1 

Pure signals can also be represented in terms of density operators. If the 
input signal \x) is deterministic (has probability one), then 



p x = \x)(x\. (4) 

Since this formalism allows to describe both types of quantum signals, in 
the following we shall characterize the signals at the input and output of our 
systems by the density operators p x and p y , respectively. As for the relation 
between the input and output signals, it can be expressed by means of a 
superoperator $: 



Py = $(Ar)- (5) 
Superoperators are complete positive operators that transform density oper- 



ators preserving its hermiticity and trace properties (see [24] for an account 
of superoperator theory). As we shall study linear quantum systems, super- 
operators will be also necessarily linear. 

In its more general form, the action of a superoperator on a quantum 
signal can be written as 

Py = J2^P X Ml, (6) 

where the operators, called Kraus operators, satisfy the closure relation 
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£MtM M = i. (7) 

It can be shown [24] that the necessary number of Kraus operators is at 
most iV 2 . Equation (|6|) is the Kraus representation of the density operator 
p y . When this representation requires only one Kraus operator, the system 
is said to be invertible since the Kraus representation reduces to a unitary 
transformation. It is remarkable, and we shall have the opportunity to give 
examples below, that any conceivable quantum system can be modelled by 
a superoperator. 



4 Matrix representation of QIP 

In the previous section we have seen that quantum signals are vectors be- 
longing to a signal space and, more generally, density operators. Up to now 
the discussion has remained rather abstract. In this section we shall develop 
a matrix formalism of the input-output theory previously presented. 

Consider an arbitrary pure signal \x) and let us define on £ the orthonor- 
mal base {\u{); i = 1, . . . , N}. In this base \x) can be written as: 

N 

\ x ) = X MK)> ( 8 ) 

n=l 

where the coefficients of the expansion are given by the inner products 
(u n \x). Given the base used, these coefficients univocally determine the 
information contained in the signal. Therefore, in a given base, the quan- 
tum signal can be represented by the sequence x[n], much as in classical 
DSP. 

Let us turn now to the more general case described by density operators. 
When the signal-space base just introduced is used on the signals \xi) of 
Section ||, we get 

N 

\Xi) = ^2 x i[ n }\ u n), (9) 
71=1 

where Xi[n] = (u n \xi). If this expansion is used in equation (|3|), the following 
representation of the input density operator is obtained: 

N 

Px = Px[n,m}\u n )(u m \, (10) 

71, 771=1 
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where 



M 

p x [n,m] =^2p i Xi[n]x*[m\. (11) 
i=l 

Now the information at the input of the system is contained in the hermitian 
matrix p x whose elements are given by equation ([ll]). From the properties 



of density operators [23] it can be demonstrated that: 

1. The eigenvalues of p x are always between and 1. 

2. The trace of p x is 1. 

3. p x is positive definite. 

Note that in the case of the pure signal \x), p x reduces to an idempotent 
matrix of elements p x [n,m] = x[n]x*[m]. 

The action of a superoperator can also be written in matrix form once a 
base of the signal space is chosen. Using the same base as above, the Kraus 
representation of the quantum output signal, see equation (||), is given by 

N 

Py [n,m] =J2J2 M M Mp^,i]Mt[j,m], (12) 

fj. i,j=i 

where M^kJ] = (u^M^ui) and, from (^), 

N 

^^Mt^ik^M = <%. (13) 

A* k=l 

Equation (|l2]) can be regarded as the equivalent of the classical input-output 
relationships (||) or (|2|) in QIP. In general, the dynamics of quantum systems 
can be described employing just standard matrix algebra: the output signal 
is obtained as the sum of all the congruence matrix transformations of the 
input signal given by the matrices associated to the Kraus operators. In 
compact form: 

A* 

Notice the more complex character of the quantum formalism as compared 
to the classical one. This complexity comes from the fact that the quantum 
description has to incorporate the possibility of quantum decoherence, i.e. 
the transformation of an input pure signal into a mixed output one. 
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In order to clarify the concepts introduced, we present a simple example. 
Consider the quantum-mechanical binary communication channel depicted 
in Figure 1. The signal space in this case is two-dimensional. If we denote 
the two elements of a base of this space by |0) and |1), then any signal 
(ordinarily called qubit, quantum bit) can be written as the superposition 
a|0) + (3\l), where \a\ 2 + |/3| 2 = 1. Let us assume that the input signal can 
be |0) or |1). According to the figure, the output signal can also be |0) or 
|1) with the probabilities shown in the figure. In this simple case the matrix 
form of one set of Kraus operators that describes the channel is: 



(15) 

These Kraus operators can also be used to obtain the output signal 
when at the input of the channel we have a signal of the more general form 
a\0)+fi\l): 



p\a 2 \ + (1 -p)\/3\ 2 pa/3* 

pa* (3 p\(3 2 \ + (1 -p)\a\ 2 



(16) 



In this case the input signal will remain the same with probability p, will 
transform to |0) with probability (1 — p)|a 2 |, and to |1) with probability 
(l-p)\(3 2 \. 



5 Composite quantum systems: Entangled signals 

Many of the proposed applications of QIP, such as Quantum Computation, 
Quantum Cryptography and Quantum Communications are based on the 
concept of information entanglement (see, e.g. pq| ). Entangling information 
in two or more signals has some advantages over the one-signal-one-system 
QIP introduced in the previous section. In this section we present com- 
posite systems. We shall show that the input-output formalism previously 
presented can also describe this kind of quantum systems. 

Consider two quantum systems A and B whose input and output signals 
are, respectively, the density operators p£, p^ , py, and p^ . Let us assume 
that the input-output relationship for each system is known (in the form of 
two superoperators $a and $#), and that there is some interaction between 
the actual physical processes modelled by the two systems. The nature of 
this interaction is not relevant to our current purpose. One can define a 
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global or composite system that comprises the two systems and its physical 
interaction, and global input and output signals. In standard QM this is 
described making use of the tensor product of signal spaces: If £a and £b 
are, respectively, the signal spaces in systems A and B, then £a ® £b is the 
signal space of the global signals. Therefore, the composite system AB has 
global input and output signals given, respectively, by p% B and py B , and its 
input-output relation is described by some global superoperator: 

P} b = $ab(p£ B )- (17) 

Like the signals and systems of the previous section, composite systems 
and global signals have also a matrix representation. If the bases of £ a and 
£b are, respectively, {|a n ); n = 1, • • • , A^} and {\b n }; n = 1, • • • , N} (we are 
assuming, for simplicity, that both signal spaces have the same dimension, 
N), the base of £ab is {\a n ) ® \b m )] n,m = 1,---,N}. In this base the 
operator density associated to the global input signal is 

N 

Px B = p£ B [ n , m ',hj]\ a n,b m }(ai,bj\. (18) 

n,m,i ,j=l 

As the number of signals and systems is doubled, so does the complexity of 
the description, and now the signals are described by rank-four tensors. No- 
tice, however, that the information contained in the tensor can be arranged 
in a standard square matrix with elements p^ B [k,l], where the indexes k 
and I (k,l = 1,---,N 2 ) enumerate, respectively, the different (A^ 2 in each 
case) pairs (n, m) and 

Although we have described the action of the composite system globally, 
in some situations it may be necessary to access the A or B part of a global 
signal. This can be accomplished using the partial trace |23f : 

N 

pj[n,i] = Tr B {pj B [n,m;i,j]} = ]T p£ B [n, m; i, m], (19) 

m=l 
N 

p B [m,j]=TrA{Py B [n,m;i,j]} = £ pf[n, m; n, j]. (20) 

n=l 

This may be the case of a physical system that interacts undesirably with 
another. Using the equations above the influence of the disturbing system 
on the dynamics of the main system can be studied. 

In some circunstances, global signals of composite systems cannot be fac- 
tored accordingly to the tensor product defined on the global signal space, 
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i.e. p AB / p A <g) p B , or, more generally, p AB / J2iPiP A ® Pf ■ This means 
that the information contained in all the partial signals is not equivalent to 
the information contained in the global signal. Thus, quantum information, 
unlike its classical counterpart, can be encoded in nonlocal correlations be- 
tween the different parts of a composite system. When this happens, signals 
are said to be entangled. 

In order to clarify the ideas presented in this section, assume that the 
signal at the input of system A is | tp x )A = a|0)/i + (3\ 1)b and that the signal 
at the input of system B is \ifi x )B = |0)_b- In the bases {10)^, and 
{|0)b) the input signals can be written, respectively, as 



Px 



| a 2 1 a/3* 
a*j3 \j3 2 \ 



and 



Px 




(21) 



(22) 



In the natural base of the composite system, {|00), |01), |10), |11)}, the global 
input signal is represented by the matrix 



pT 



( \a 2 \ af5* 



(3a* \P\ 2 

\ 






0/ 



(23) 



Consider now that the composite system AB produces the unitary tranfor- 
mation p AB = M/9^ B M^ on the global input signal, where 



M 



(24) 



/ 1 \ 
10 
1 
\ 1 J 

For the input given by equation (^), the corresponding global output signal 
is 



Pf 





\a\ 2 








a(3* \ 






























V 


a*P 








l/?| 2 ) 



(25) 
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It can be shown that, contrary to what happened at the input, at the output 
the signal cannot be expressed as the tensor product of a signal from Ea and 
another from Eb- The signals, that were uncorrelated at the input, are now 
entangled. 



6 Some quantum systems 

In Section || we have shown that the operation of any quantum system can 
be described in matrix form (see equation (14)). In this section we are going 
to show that the matrix formalism proposed can in fact be applied to the 
description of two typically quantum processes: quantum evolution, and 
quantum measurement. 



6.1 Time evolution described as a system operation 

Physically, any kind of processing requires a lapse of time. In the systems 
previously considered we have assumed that processes occur instantaneously, 
so no time dependence has been included in the formalism. Consider now 
the case of a quantum system whose input is given by p x and its output, 
the processed input, by p y . If the processing begins at the instant to and 
finishes at t, the laws of QM say (see, e.g. p6| ) that both quantum signals 
are related by the unitary transformation 

p y = \J(Ar)p x V^Ar), (26) 

where U is the matrix representation of the standard evolution operator 
U = exp^—iHAr/fi), with H the Hamiltonian of the system, and At = 
t — to- Equation (pq), owing to the unitarity of the matrix U, is a particular 
case of equation fll4j) , and thus can be seen as a Kraus representation. 

The evolution represented by equation (^6|) is characteristic of isolated 
quantum processes, i.e. closed physical systems. However, in some cases 
two or more processes are coupled forming a composite physical system and 
one is interested in the time evolution of the global signal and its different 
parts. This composite physical system can be regarded as a closed system, 
but the parts (the so-called open systems) certainly not. This situation, as 
we show below, can also be described by the matrix formalism. 

As in Section ||, assume we have two interacting physical systems, A and 
B, that unite to form the composite physical system AB. This system is 
closed, so it can be described by an equation like p^): 
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pf = \J AB (Ar)p^ AB (Ar), (27) 

where the Hamiltonian inside U takes into account the dynamics of both 
physical systems as well as their mutual interaction. Taking the partial 
trace, the matrix representation of the time-evolution of one of the partial 
systems, say A, can be easily obtained: Py = Ttsipy 3 }- However, it 
is slightly more involved to obtain the Kraus representation of this input- 
output relationship. If the two systems are initially uncorrelated, i.e. = 
Px®Pxi it can be shown [27] that 



N 

p* = J2 "»pM> (28) 

71=1 



where the fi n matrices, that obey the condition (13), are function of the 



evolution of the composite system, Uab(At), and of p^ 



6.2 Measurement described as a system operation 

The quantum description of measurement is still rather controversial, par- 
ticularly among researchers focused on the interpretative aspects of QM 
[p8| . Measurement can be described ideally as a process that takes place 
instantaneously (the von Neumann or othogonal measurement) in a closed 
system, or taking into account that any measurement needs an interaction, 
and that interactions do not occur (at least in the physical world) instan- 
taneously (the generalized measurement in an open system) p9| . We shall 
show that the general formalism developed in Section |I| can describe both 
kinds of measurement processes. 

Consider a quantum system that performs a measurement on its input 
signal p x to produce at its output, as a result of such a measurement, the 
signal p y . Associated to any measurement there is an observable Q (an her- 
mitian operator; see, e.g. |^]) with the eigenvalue problem (for notational 
simplicity we study the non-degenerate case) 

Q\ln) = q n \q n )- (29) 
The eigenstates of this problem are a base of the signal space, and the family 
of orthogonal projectors {\q n ) (<?n| ; n = 1, • • • , iV} are such that J2n=i \ln) (<?n| 
1. It can be shown from QM p6| that, in operator language, 

N 

Py = ^(l^>(^l)Ar(l<7n>(9nl)- (30) 

ra=l 
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This equation can be written in matrix form once a base of the signal space 
is selected: 



N 

Py = Y, M nP x Ml (31) 
n=l 

where the elements of the Kraus matrices are given by M n = (ai\q n )(q n \ aj ) 
and we have used the base {\a n ),n = 1, • • • , N}. 

As we mentioned at the beginning of this subsection, any physical mea- 
surement requires the interaction of the signal with (at least) the measuring 
apparatus, and this interaction is not an instantaneous physical process. 
The quantum-mechanical consideration of this fact leads to a more complex 
description of measurement based on Positive Operator Valued Measures 



(POVM) [3C]. In these generalized measurements the signal to be mea- 
sured interacts with another called ancilla and an orthogonal measurement 
is performed on the ancilla part of the resulting global signal. Owing to 
the interrelation between both signals, this measurement gives information 
about the original signal. 

Suppose we have a signal p^ that interacts with another signal another 
initially uncorrelated signal p B . The time evolution of this interaction, ac- 
cording to what we saw in Subsection |6.1| , can be modelled by the unitary 
Kraus representation (]27|). Let us now perform the orthogonal measure- 
ment of Q on the B part of the output signal Py B , and leave the A part 
untouched. According to the previous formalism for orthogonal measure- 
ments, see equation (|30[), the resulting composite signal is 

Pz B = E(|«fc><?*| ® i A )UAB(Ar)^ B U\ B (Ar)(\ qk )( qk \ l A ). (32) 

k=l 

If we now take the partial trace over B, the output signal of system A after 
the generalized measurement can be written as 

pf = Tr B {U AB p£ B U\ B }. (33) 

An input-output relationship of the form of equation (^) can be obtained 
for the A system, just as we did in Section |6.1| . 

7 Conclusion 

We have shown that any quantum-mechanical information processing can 
be seen as the operation of an abstract system with input and output sig- 
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nals. Both signals and systems are represented by matrices. The output 
signal is obtained as a superposition of congruence matrix transformations 
on the input signal. Our main objective in proposing such a formalism is to 
provide a framework to further study QIP suitable for non-physicists. As 
two applications of the formalism, we have expressed in signal-system form 
the processes of quantum time-evolution and quantum measurement. 
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